Coupling molecular spin states by photon-assisted tunneling 
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Artificial molecules containing just one or two electrons provide a powerful platform for studies of 
orbital and spin quantum dynamics in nanoscale devices. A well-known example of these dynamics 
is tunneling of electrons between two coupled quantum dots triggered by microwave irradiation. So 
far, these tunneling processes have been treated as electric dipole-allowed spin-conserving events. 
Here we report that microwaves can also excite tunneling transitions between states with different 
spin. In this work, the dominant mechanism responsible for violation of spin conservation is the spin- 
orbit interaction. These transitions make it possible to perform detailed microwave spectroscopy of 
the molecular spin states of an artificial hydrogen molecule and open up the possibility of realizing 
full quantum control of a two spin system via microwave excitation. 



In recent years, artificial molecules in mesoscopic sys- 
tems have drawn much attention due to a fundamental in- 
terest in their quantum properties and their potential for 
quantum information applications. Arguably, the most 
flexible and tunable artificial molecule consists of cou- 
pled semiconductor quantum dots that are denned in a 
2-dimensional electron gas using a set of patterned elec- 
trostatic depletion gates. Electron spins in such quan- 
tum dots exhibit coherence times up to 200 fis [1], about 
10 4 — 10 6 times longer than the relevant quantum gate 
operations 0], making them attractive quantum bit 
(qubit) systems [j. 

The molecular orbital structure of these artificial quan- 
tum objects can be probed spectroscopically by mi- 
crowave modulation of the voltage applied to one of the 
gates that define the dots [5[. In this way, the delocal- 
ized nature of the electronic eigenstates of an artificial 
hydrogen-like molecule was observed [H, 0- More re- 
cently, electrical microwave excitation was used for spec- 
troscopy of single spins [Mil and coherent single-spin 
control [a,llflj, via electric dipole spin resonance (EDSR). 

Here, we perform microwave spectroscopy 0,0, EH on 
molecular spin states in an artificial hydrogen molecule 
formed by a double quantum dot (DD) which contains 
exactly two electrons. In contrast to all previous PAT ex- 
periments, we observe not only the usual spin-conserving 
tunnel transitions, but also transitions between molecu- 
lar states with different spin quantum numbers. We dis- 
cuss several possible mechanisms and conclude from our 
analysis that these transitions become allowed predomi- 
nantly through spin-orbit (SO) interaction. The possibil- 
ity to excite spin-flip tunneling transitions lifts existing 
restrictions in our thinking about quantum control and 
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detection of spins in quantum dots, and allows universal 
control of spin qubit s without gate voltage pulses. 



1. DEVICE AND EXCITATION PROTOCOL 

Fig. la displays a scanning electron micrograph of 
a sample similar to that used in the experiments. It 
shows the metal gate pattern that electrostatically de- 
fines a DD and a quantum point contact (QPC) within 
a GaAs/(Al,Ga)As two-dimensional electron gas. An 
on-chip Co micro-magnet (/imagnet) indicated in blue 
in Fig. la generates an inhomogeneous magnetic field 
across the DD, which adds to the homogeneous external 
in-plane magnetic field B (see the appendix A for more 
sample details), but is not needed for the molecular spin 
spectroscopy. The sample was mounted in a dilution re- 
frigerator equipped with high-frequency lines. The gate 
voltages are set so that the DD can be considered as a 
closed system (the inter dot tunneling rates are 10 4 times 
larger than the dot-to- lead tunneling rates), and the tilt 
of the DD potential is tuned by the dc-voltages Vl and 
Vr, applied to the left and right side gates. Working near 
the turn-on of the first conductance plateau, the current 
through the QPC, Iqpc, depends upon the local charge 
configuration and provides a sensitive meter for the ab- 
solute number of electrons (n^^i?) in the left and right 
dot, respectively [l2l [l3|. 

First, we excite the DD as indicated in the top panel 
of Fig. lb, by adding to Vr continuous- wave microwave 
excitation at fixed frequency v = 11 GHz. When the 
photon energy of the microwaves matches the energy 
splitting between the ground state and a state with a 
different charge configuration, a new steady-state charge 
configuration results, which is visible as a change in the 
QPC current, AIqpc- The excitation is on-off modu- 
lated at 880 Hz and lock- in detection of AIqpc reveals 
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FIG. 1: Photon-assisted tunneling in a 2-electron dou- 
ble quantum dot. a, Scanning-electron micrograph top 
view of the double dot gate structure with Co micromagnet 
(blue). The voltages applied to the left Vl and right Vr side 
gates (red) control the detuning e of the double dot potential. 
The double dot charge state is read out by means of the cur- 
rent Iqpc running through a nearby quantum point contact 
(white arrow), b, Charge stability diagram around the 2- 
electron regime at B = 1.5 T. (til,tir) indicate the absolute 
numbers of electrons in the left and right dot, respectively. 
During measurements, 11 GHz microwaves with 880 Hz on- 
off modulation are applied to the right side gate. The top 
panel displays schematically one cycle of the ac signal applied 
to the right side gate. Along the detuning axis (dashed ar- 
row) PAT-lines are observed, c, In the conventional picture 
of PAT, the first sidebands seen in b should appear when the 
detuning of the (0,2) and (1,1) states matches the photon 
energy and interdot tunneling is induced. Further sidebands 
are then interpreted as multi-photon transitions. fiF is the 
chemical potential of the left and right electron reservoir, d, 
Same as in b, but the microwaves are interrupted every 5 fis 
by a 200 ns, P £ — 2 mV detuning pulse applied to the left 
and right side gates (see the schematic in the top panel). The 
pulses generate a reference line (black arrow) due to mixing 
at the ST + anti-crossing. 



the microwave- induced change of the charge configura- 
tion (see the appendix A for further experimental de- 
tails). The lower panel of Fig. lb shows AIqpc as a 
function of Vl and Vr near the (1,1) to (0, 2) boundary 
of the charge stability diagram. Sharp red (blue) lines in- 
dicate microwave-induced tunneling of an electron from 
the right to the left dot (left to right), labeled as An = +1 
(An = — 1), respectively (see Fig. lc). Sidebands can re- 
sult from multi-photon absorption. At the boundaries 
with the (0, 1) and (1,2) charge states, no energy quan- 
tization is observed, since here electrons tunnel to and 



from the electron- state continuum of the leads. At first 
sight, the observations in Fig. lb thus appear to be well 
explained by the usual spin-conserving PAT processes. 

Surprisingly, the position in gate voltage of the reso- 
nant lines exhibits a striking dependence on the in-plane 
magnetic field, B. This is clearly seen in Figs. 2a and 2b, 
which display the measured PAT spectrum along the DD 
detuning e axis (dashed black arrow in Fig. lb) as a func- 
tion of B for 20 GHz and 11 GHz excitation, respectively. 
Since the gate constitutes an open-ended termination of 
the transmission line, the excitation produces negligible 
AC magnetic fields at the DD, and is therefore expected 
to give rise to only electric-dipole allowed spin-conserving 
transitions, with no B dependence. Furthermore, there 
is a pronounced asymmetry between the position of the 
red and blue PAT lines. 

In these figures, the detuning axis was calibrated for 
all magnetic fields by introducing a reference line (see 
also Fig. Id, black arrow) that facilitates interpretation 
of the spectra despite residual orbital effects of the mag- 
netic field. This line was produced by interspersing the 
microwaves every 5 fis with 200 ns gate voltage pulses 
along the detuning axis (see top panel in Fig. Id), lead- 
ing to singlet-triplet mixing as described in Ref. [2] . The 
short gate voltage pulses do not noticeably alter the po- 
sition of the PAT lines (compare Figs, lb and Id). The 
reference peaks visible at around e = 200 jaeV in Figs. 2a 
and 2b were aligned by shifting all data points at a given 
B by the same amount in detuning (see the appendix for 
the full details of this post-processing step). 



2. INTERPRETATION OF THE 
PHOTON-ASSISTED TUNNELING SPECTRA 

The complexity of the PAT spectra shown in Fig. 2a, b 
can be understood in detail if we allow for non-spin con- 
serving transitions. The two diagrams in Fig. 2c show 
the energies of all relevant DD-states (four (1, l)-states 
and one singlet 5(0, 2)-state) as a function of e for two 
different (fixed) magnetic fields, i.e. the spectrum of the 
DD along the two horizontal dotted lines in Fig. 2b [lij . 
Note that the only difference between the two diagrams 
is the splitting between the three triplet T(l, l)-states. 

First we explain the resonances observed along the up- 
per dotted line in Fig. 2b (B = 2.5 T). In the cor- 
responding (upper) diagram in Fig. 2c, we plot the 
ground state energy for all e with a thick line. If the mi- 
crowave excitation is off-resonance with all transitions, 
the system will be in this ground state. For instance, at 
e = 150 /ieV, there is no state available 11 GHz above the 
5(0, 2) ground state (gray arrow) and the system stays in 
5(0, 2). However, when decreasing the detuning, at some 
point T + (l,l) becomes energetically accessible (red ar- 
row) and, since we allow for non-spin-conserving transi- 
tions, is populated due to the microwave excitation. For 
this PAT transition, the spin projection on the quan- 
tization axis is changed by Am = +1. The resulting 
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FIG. 2: Photon-assisted-tunneling spectra and simulations. a,b, Microwave induced change of the QPC current 
AIqpc as a function of the double dot detuning e and the external magnetic field B for 20 GHz (a) and 11 GHz (b) frequency, 
respectively, singlet-triplet mixing due to 2 mV detuning pulses generates a reference signal that is used to calibrate the 
detuning axis (see lower panel in Fig. 2f). c, Eigenenergies vs. double dot detuning e of the 2-electron spin states in the (1, 1) 
and (0, 2) charge regime for two external magnetic fields B — 2.5 T (upper panel) and B — 1.5 T (lower panel), respectively. 
5(0, 2), 5(1, 1) and T(l, 1) character of the eigenstates is indicated by blue, green and red color, respectively. The molecular 
spin ground state is indicated by thick lines. The vertical arrows indicate PAT transitions for a constant microwave frequency 
involving spin flips. The transition indicated by a dashed arrow is suppressed, because the initial state lies above the ground 
state. The red circle in the lower panel indicates the detuning position of the reference signal, that is generated by a detuning 
pulse with amplitude P £ to the 5T + ant i- crossing. d,e, Simulated PAT spectra for 20 GHz (a) and 11 GHz (b) frequency, 
respectively. The color indicates the change of the population of the steady-state charge state An as would be observed in 
an on-off lock-in detection. A finite temperature of 100 mK and spontaneous relaxation via the phonon bath are taken into 
account, f, Iqpc scanned with higher resolution in the anti-crossing region (black rectangle in Fig. 2b) at 11 GHz excitation. 
The green dashed lines indicate the expected detuning positions of the PAT transitions. The horizontal black dashed line 
indicates the magnetic field, at which the electron spin resonance condition is fulfilled E z = gjj>B(B + bo) = hv. The graph is 
concatenated from two scans that overlap at 1.9 T. The inset displays three magnetic fields (blue dots), at which the center of 
the horizontal blue triplet resonance line is observed, as a function of the microwave frequency v. The dashed black line gives 
the expected position of the triplet resonance. 



change of steady-state charge population (increased pop- 
ulation of (1,1), or AIqpc > 0) is detected by the QPC 
and yields the red peak in Fig. 2b. Decreasing e fur- 
ther, there are two more resonances detectable: (i) the 
5(0, 2)-5(l, 1) transition (dotted red arrow, Am = 0), 
although the signal will be weakened due to the fact that 
5(0, 2) is not unambiguously the ground state anymore. 
Note that a transition 5(0, 2)-T°(l, 1) could appear in 
nearly the same detuning position, as will be discussed 
below, (ii) the T + (l, l)-5 transition (blue arrow), where 
5 stands for the hybridized 5(0,2) — 5(1, 1) singlet, re- 
sults in a negative (blue, AIqpc < 0,Am = —1) signal 
from the charge detector since the ground state is now 
(1,1) and the excited state is (0,2). We see that this 
simple analysis explains both the positions and the signs 
of the resonances observed in the data. 

A similar analysis can be made for other magnetic 



fields. For instance, for the spectrum plotted in the 
lower diagram of Fig. 2c we find two resonances with 
A/qpc > (red arrows), and one with AIqpc < (blue 
arrow). Note that the 'blue' transition now connects the 
ground state to the other branch of the hybridized 5 com- 
pared to the high magnetic field case. Indeed, the singlet 
anti-crossing is directly probed, resulting in the two blue 
curved lines observed in the data around e — (Fig. 
2b). The fading out of the blue signal at low fields can 
be understood from pumping into the metastable state 
5(1, 1): the microwaves excite the system from T + (l, 1) 
to 5(0, 2), from where it relaxes quickly to 5(1, 1). How- 
ever, relaxation from 5(1, 1) back to the ground state is 
slow due to the small energy difference of this transition 
and the small phonon density of states at low energies 
[3 [HI- This pumping weakens the detector signal, since 
5(1, 1) has the same charge configuration as the ground 
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state. 

In order to verify this interpretation, we calculate in 
Figs. 2d,e the position and intensity of the spectral 
lines at fixed microwave frequency, based on the en- 
ergy level diagram of Fig. 2c. In the simulations, all 
single-photon transitions between the ground state and 
the excited states are allowed by including a matrix el- 
ement |T ± (1, 1))(6'(0, 2)| (see appendix A). The input 
parameters for the calculation of the resonant positions 
are the interdot tunnel coupling t c , the absolute elec- 
tron ^-factor \g\, a magnetic-field contribution bo from 
the /imagnet parallel to B as well as a magnetic field 
gradient AB = (A£^, 0, AB^) between the dots (in Fig. 
3, we show how t c , \g\ and bo can be extracted from 
the experimental spectra) . The color scale represents the 
calculated steady-state An that results from microwave 
excitation, orbital hybridization and phonon absorption 
and emission at 100 mK. All the PAT transitions visible 
in the simulation also appear in the experiment, with ex- 
cellent agreement in both the position and relative inten- 
sity of the spectral lines. Especially, the vanishing signal 
due to spin pumping is also predicted by the simulations 
which include phonon relaxation. 

When we zoom in on the boxed region of Fig. 2b, 
we see an additional horizontal blue feature at B « 2 T 
(Fig. 2f) that also appears in the calculated spectra of 
Fig. 2e. This feature is due to a triplet resonance from 
T + (l,l) to T°(l, 1) that becomes detectable by relax- 
ation into the meta-stable 5(0, 2) state. In the detuning 
range where this line appears, the 5(0, 2) state lies en- 
ergetically only slightly above the T + (l,l) state, so re- 
laxation back to the T + (l, 1) ground state is suppressed, 
again by the small phonon density of states at low ener- 
gies (see Fig. S2a of the appendix). The triplet resonance 
is expected to appear at E z = qiib(B + &o) = his, where 
h is Planck's constant and \±b the Bohr magneton. The 
inset of Fig. 2f shows the magnetic fields correspond- 
ing to the center of the measured triplet resonance line 
for three excitation frequencies (see Figs. S2b,c of the 
appendix for the spectra), which are in good agreement 
with the expected positions (black dashed lines in Fig. 
2f) based on the values \g\ and bo determined in the next 
section from other features of the PAT spectra. Sur- 
prisingly, the measured triplet resonance exhibits a finite 
slope in the B(e) spectra. A longitudinal magnetic field 
gradient AB^ gives rise to such a detuning dependence, 
but the AB\\ required in our simulations to reproduce 
the observed slope is AB^ > 80 mT/50 nm, an order of 
magnitude larger than the gradient we calculate for the 
/imagnet. The magnitude of the slope remains a puzzle. 



3. EXTRACTING ARTIFICIAL MOLECULE 
PARAMETERS 

We now show how |g|, t c and 6q? the parameters used 
for all simulations, can be extracted independently from 
the experimental spin-flip PAT spectra. For this analysis 
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FIG. 3: Analysis of the photon-assisted-tunneling 
spectra, a, The detuning difference As between the spin 
conserving line Am = and the Am = dzl line is plotted as 
a function of the external magnetic field B at 20 GHz (see 
Fig 2a), in order to fit the absolute effective electron g- factor 
g from the slopes of the linear fits (solid lines), inset, Ae for 
the blue Am = +1 line at 11 GHz for different tunnel cou- 
plings. The offset of the curves increases with increasing t c 
(red to black circles) while g remains constant, b, As between 
the Am = line and the Am = —1 (1,1) to (0, 2) transition 
from Fig. 2c. The fit of the anti-crossing (red line) allows for 
a precise determination of t c . 



we only use the relative distance Ae between PAT lines 
at fixed magnetic field, in order to be independent from 
the calibration of the detuning axis by means of the ref- 
erence line. Fig. 3a shows Ae ± as a function of B using 
the 20 GHz data. Ae ± is defined as the difference in 
detuning between the red Am = ±1 and Am = PAT 
lines. For a fixed both Ae ± increase linearly with the 
Zeeman energy and therefore allow fitting of \g\. (Note 
that for Ae + the linearity is only exact for sufficiently 
large B, at which the singlet anti-crossing does not af- 
fect the T + (l, 1) energy; see the appendix for a detailed 
discussion). A least-squares fit to the Ae + data gives 
\g\ = 0.382 ± 0.004 (Fig. 3a). From the linear behavior 
of Ae + , we also deduce that there is negligible dynamic 
nuclear polarization in the experiment. 

Knowing \g\ precisely, we make use of the blue anti- 
crossing in Figs. 2b and 2f in order to determine t c (and 
bo). Fig. 3b shows the difference in detuning Ae' between 
the blue Am = — 1 and the red Am = lines in Fig. 2f. 
Assuming the Am = line corresponds to the 5(0, 2)- 
5(1, 1) transition (as shown below), then 
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where the first term is the detuning position of the 
T + (l,l)-5 transition and the second the one of the 
5(0, 2)-5(l, 1) transition. The best fits are obtained with 
t c = 8.7 ± 0.1 /xeV and b = 109 ± 16 mT. This value for 
bo matches very well our simulations of the stray field of 
the /imagnet at the DD location (see appendix A). 

An important question left open so far is whether 
the red Am = line involves predominantly transitions 
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from 5(0,2) to 5(1,1) or to T°(l,l). The transition 
to 5(1,1) does not require a change in the (total) spin 
and is thus expected to be excited more strongly than 
that to T°(l,l). However, relaxation from 5(1,1) back 
to 5(0, 2) will be stronger as well, so it is not obvious 
what steady-state populations will result in either case. 
Furthermore, given the small energy difference between 
5(1, 1) and T°(l, 1), the two transitions are not resolved 
in Fig. 2. Fig. 3a helps to answer this question: The 
observation that As+ > As~ indicates that the Am = 
line originates from the transition to 5(1,1) and not to 
T°(l, 1). For the former we expect Ae ± = E z ± J, with 
t 2 

J(e,t c ) = -f + 0(tc) the exchange energy, whereas the 
latter would result in Ae + < Ae~ (in both scenario's, bo 
causes an additional fixed offset in both Ae ± , but it does 
not contribute to their difference). This interpretation is 
consistent with the increase of Ae + with larger interdot 
tunnel coupling, hence larger J (Fig. 3a inset; note that 
the slopes are not affected). It is further supported by 
the data in Fig. S3b. 

So far only single-photon processes were considered, 
but at higher microwave power, also multi-photon lines 
emerge (Fig. 4a), mostly for the 5(0, 2)-5(l, 1) transi- 
tion (green dashed lines in Fig. 4a). Like the single- 
photon 5(0, 2)-5(l, 1) line, their position in detuning is 
^-independent (see appendix for details). 



4. IDENTIFICATION OF THE SPIN-FLIP 
MECHANISMS 

Having shown the power of spin-flip PAT for detailed 
molecular spin spectroscopy, we now discuss the mech- 
anisms responsible for this process as confirmed by our 
simulations. As a first possibility, the transitions from 
5(0, 2) to the triplet (1, 1) states can take place through a 
virtual process involving 5(1, 1): the state 5(0, 2) is cou- 
pled to 5(1, 1) by the interdot tunnel coupling, and an 
(effective) magnetic field gradient AB = (AB^, 0, AB^) 
across the DD couples the spin part of all the (1,1) states 
to each other [l3|, [l6|. Here AB has a contribution 
from the effective nuclear field and from the /imagnet. 
The transition matrix element from 5(0, 2) to T ± (l, 1) is 

oc t c A ^ , assuming E z ^> J. In the following, we use the 
independence of this process as a fingerprint and focus 
on the red Am = +1 line in Figs. 2a and 2b, as we can 
follow it over the entire magnetic field range. The inten- 
sity of this line is constant in B and even if the microwave 
amplitude E is varied, we observe no ^-dependence in 
the area under this peak (Fig. 4b). Before we conclude 
that the transition rate is magnetic field independent, we 
recall that the observed PAT lines reflect the steady-state 
change in the charge configuration resulting from stim- 
ulated photon emission and absorption and spontaneous 
relaxation. In order to rule out that a field-independent 
steady state is reached from a field dependence of re- 
laxation and excitation that cancel each other, we verify 
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FIG. 4: Power-dependence of the photon-assisted tun- 
neling spectra and spontaneous relaxation, a AIqpc 
as a function of the double dot detuning e and the microwave 
amplitude E for 11 GHz and B = 1.0 T (middle panel) and 
B = 2.5 T (lower panel). The excitation scheme from Fig. 
Id is employed with reference pulse amplitude P £ — 1.5 mV. 
The green (red) dashed lines mark the multi-photon Am = 
(Am = +1) PAT transitions. The reference signal stemming 
from the pulse is marked by orange arrows. The voltage am- 
plitude E is measured at the end of the coaxial lines at room 
temperature. The uppermost panel displays a linecut mea- 
sured at 1 T. The red line is a least-squares fit by a sum of 
four lorentzian peaks to the data, b The fitted Lorentzian 
peak area of the Am = +1 line from a is plotted as a func- 
tion of microwave amplitude E and magnetic field B. c The 
fitted Lorentzian linewidth at half maximum (FWHM) of the 
Am = +1 and the lines Am = from a is plotted as a func- 
tion of microwave amplitude E for B — 1 T. The error bars 
in b,c and are determined from the Lorentzian least-squares 
fit (see the inset of Fig. 4a). 



that the spontaneous relaxation rate is field-independent 
as well (see appendix). These observations suggest that 
the coupling mechanism is magnetic field independent 
and thus virtual processes involving 5(1, 1) do not give a 
strong contribution to the transition rates. 

More recently, two mechanisms were considered that 
provide a direct, ^-independent matrix element between 
5(0, 2) and the (1, 1) triplet states: (i) the hyperfine con- 
tact Hamiltonian is of the form . S(r — rj)P • 5. Thus, 

nuclear spins, Ij, in the barrier regions, where the 5(0, 2) 
spatially overlaps with each of the (1,1) triplet states, can 
flip-flop with the electron spin, 5, simultaneously with 
charge tunneling [17]. (ii) The SO Hamiltonian is of the 
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form p X) yS x ^ y and can directly couple states which differ 
in both orbital and spin 

BO- (When the orbital part 
of the initial and final state are the same, the SO Hamil- 
tonian does not provide a direct matrix element and the 
transition rate becomes Independent [8|, |2Ch23| .) The 
ratio of the SO mediated rate and the hyperfine medi- 
ated rate can be estimated as ^tt^ t^~ ? which is a few 

4A Aso ' 

thousand in the experiment (see the appendix). Here 
Eq is the single-dot level spacing, N the number of nu- 
clei in contact with one dot, A the hyperfine coupling 
strength, d the interdot distance and Xso the spin-orbit 
length. We therefore believe that SO interaction is the 
dominant spin- flip mechanism for the observed PAT tran- 
sitions. The presence of a magnetic field independent 
matrix element between 5(0,2) and the (1, 1) triplets is 
confirmed by the observation that the intensity of the 
reference signal shows no field dependence. 

Finally, we extract from Fig. 4a the fitted linewidth as 
a function of driving power (Fig. 4c). For small E, we 
find a width > 1 GHz, similar to that observed earlier for 
spin-conserving PAT processes [6]. For stronger driving, 
both the Am = and the Am = +1 lines are further 
broadened, up to E « 1.5 mV. If these lines were power 
broadened, their width would imply transition rates in 
excess of 1 GHz. However, we do not believe that this 
is the case, since in measurements with short microwave 
bursts the populations saturated only on a long (10 [is) 
timescale (data not shown). Presumably charge or gate 
voltage noise is responsible for the broadening instead. 

We have shown that in our DD system, all (1,1) spin 
states have at least weakly allowed electric dipole tran- 
sitions to 5(0,2). In materials with high SO interaction 
like In As, the effect of the non-spin conserving PAT will 
be even stronger. In materials with weak SO interac- 
tion, a strong gradient magnetic field can be used to fa- 
cilitate spin-flip tunneling transitions. In all cases, this 
opens the possibility of realizing full quantum control of 
the (1,1) spin space via off-resonant (microwave) Raman 
transitions through the excited 5(0, 2) state, which en- 
ables a variety of new approaches to manipulating and 
even defining qubits in DDs. Furthermore, such control 
enables new measurement techniques that do not rely on 
Pauli spin blockade [24|. An example is a measurement 
that distinguishes parallel from anti-parallel spins while 
acting non-destructively on the 5(1, 1) — T°(l,l) sub- 
space, by coupling resonantly the T + (l, 1) and T~(l, 1) 
states to 5(0, 2) followed by charge readout. This con- 
stitutes a partial Bell measurement and leads to a new 
method for producing and purifying entangled spin states 

Appendix A: Methods 

1. Sample fabrication 

30 nm thick TiAu gates are fabricated on a 90 nm 
deep (Alo.3,Gao.7)As/GaAs two-dimensional electron gas 



(2DEG) by means of ebeam lithography. The double 
dot axis is aligned along the [110] GaAs crystal direc- 
tion (z-direction), which is parallel to the external mag- 
netic field direction B (Fig. la). The 2DEG is Si 5- 
doped (40 nm away from the hetero-interface), exhibits 
an electron density of 2.05 x 10 11 cm -2 and a mobil- 
ity of 2.06 x 10 6 cm 2 /Vs at 1 K in the dark. The 
grounded, 275 nm thick, 2 fim wide and 10 fim long Co 
/imagnet is evaporated on top of a 80 nm thick dielectric 
layer, aligned along B (magnetic easy axis) and placed 
« 400 nm away from the closest dot center. We calcu- 
late (26| that at the double dot position the /imagnet adds 
bo « 110 mT to B z and generates a magnetic field gradi- 
ent of AB\\ ~ 6 mT/50 nm and a transverse gradient of 
AB 1 ^ « -6 mT/50 nm at saturation (B z > 2 T). 



2. Measurement 

The sample is mounted in an Oxford KelvinOx 300 
dilution refrigerator at 30 mK. Left and right side gate 
voltages, Vl and Vr, are set by low-pass filtered dc lines 
and « 60 dB attenuated coaxial lines combined with bias- 
tees with a cutoff frequency of 30 Hz. The pre-amplified 
current through the quantum-point contact is read out 
by a lock-in amplifier locked to the 880 Hz on-off modula- 
tion of the microwaves. The bias across the double dot is 
set to /iV. Voltage pulses to the left and right side gates 
are generated with a Sony Textronix AWG520. The mi- 
crowaves are generated with a HP83650A and combined 
with the pulses to the right side gate. Microwave bursts 
and detuning pulses are synchronized to ensure that the 
microwave excitation is switched off during the detuning 
pulses that generate the reference signal (see Fig. Id). 



3. Simulation 

The Hamiltonian describing the two-spin system near 
the (1, l)-(0, 2) transition is taken to be a five-state sys- 
tem, with four (1, 1) spin states and a (0,2) spin singlet 
[l6j in the presence of an external magnetic field B and 
a magnetic field gradient A£?, which includes both the 
quasi-static nuclear field and the field from the /imagnet. 
This is given by H = gfJ, B Pii(Si • (B + AB) + 5 2 • (B - 
AB))P n -e\S(0, 2))(5(0, 2)| +H U were P n is the projec- 
tor onto the (1,1) subspace, e is the detuning due to the 
difference in gate potentials from the left and right gates 
and the tunnel coupling H t = t c (\S(l, 1)) (5(0, 2)|) + 
*so(|T+(l, 1)) (5(0, 2)|) + t SO (\T-(l, 1)) (5(0, 2)|) + h.c. 
with t c the spin-conserving tunnel coupling and tso the 
spin-orbit coupling set to 5 % of t c . 

To find the signal we expect theoretically from the ex- 
periment, we add a weak, rapidly oscillating term to the 
Hamiltonian: e — >• So + Qcos(vt). We diagonalize H 
with Q = 0, then make a rotating frame transformation 
in which levels are grouped into bands n (defined by a 
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projector P n ) where the states in a band n are much 
closer in energy than his, while the energy difference be- 
tween states in band n and n + 1 are within 2/3rds of 
hv. Each band rotates at a rate nv, and we can then 
make a rotating wave approximation, keeping terms due 
to 5 that couple adjacent bands, i.e., our perturbation in 
the rotating frame and rotating wave approximation is 
V = ^En P n|(0,2)^)((0,2)5|P n+1 +h.c. Next, we add 
dissipation and dephasing by including relaxation due to 
coupling of the electron charge to piezoelectric phonons 
in a two-orbital (Heitler-London-like) model. 



Appendix B: Calibration of the detuning axis 

The photon-assisted tunneling (PAT) spectra in Figs. 
2a and 2b of the main article are measured both with a 
high energy resolution along the double dot (DD) detun- 
ing axis and over a wide external magnetic field range. In 
the experiment, we observe a monotonous, reproducible 
drift of the stable charge regions predominantly along 
the right side gate voltage as we change the magnetic 
field. Changing the voltages applied to the left and to 
the right side gate accordingly, we partially compensate 
for this drift. We then record an 11 GHz PAT spectrum 
as displayed in Fig. SI a. In order to precisely calibrate 
the detuning axis for all B, a reference signal is generated 
together with the PAT spectrum by interspersing the mi- 
crowaves every 5 /is by a 200 ns detuning pulse with am- 
plitude P £ towards negative detuning (Fig. Sib). 200 ns 
are found to be sufficient to mix the 5(0, 2) state entirely 
with the T + (l,l) state at their anti-crossing £st+, so 
that a Pauli spin blocked T + (l, 1) signal is observed at a 
detuning e = £st+ + Pe- The magnitude of P £ is chosen 
such that the reference signal appears at a detuning po- 
sition far away from the PAT signal. The pulses do not 
alter the detuning position of the PAT resonances. In 
addition to ST + mixing, mixing of the S with T°(l, 1) 
and T~(l, 1) is observed due to the pulsing. The former 
gives rise to a positive AIqpc background on the left of 
the reference signal in Fig. Sib. The latter generates a 
weak second reference line that overlaps with the 5T + 
reference line for B — ^ T, but shifts towards negative 
detuning as B is increased. 

In a post-processing step, we separately fit the position 
of the 5T + peaks for all B and shift every row of the 
spectrum, such that the peaks are vertically aligned at 
£* = P £ as shown in Fig. Sic. Thus, e* = is at the 
5T + anti-crossing, by definition. All data points at a 
given B are shifted by the same amount in detuning. 
This £*-detuning axis is well-defined but 'moves' with 
respect to the ^-detuning axis as a function of B, since 
e* = e — £$T+(B). The lever arm for the voltage to 
energy conversion is read from the voltage distance of 
the second and third Am = PAT line at v = 11 GHz 
(see Fig. 4a), which equals the photon energy hv [7]. For 
the analysis of all spectra we used the same lever arm. 

For better readability of the PAT spectra, we finally 
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Fig. S 1: Calibration of the detuning axis for the PAT 
spectra, a, PAT spectrum as measured without pulses (see 
inset of Fig. lb) with a microwave frequency of 11 GHz. The 
voltage applied to the right side gate Vr and left side gate Vl 
is swept simultaneously to follow the detuning axis indicated 
in Fig lb. b, same PAT spectrum, but with 200 ns pulses that 
intersperse the microwave bursts every 5 /is (see top panel 
of Fig. Id). A red reference line appears at positive DD 
detuning, when the detuning pulse reaches exactly the ST + 
anti-crossing. Also a weak line due to mixing at the ST~ 
anti-crossing is observed. Both lines are marked by arrows. 
The more red background is due to ST° mixing, c, PAT 
spectrum with calibrated detuning axis. The horizontal lines 
of the spectrum are shifted so that the reference pulse appears 
at the detuning position, which equals the pulse amplitude P £ 
converted into energy (see text) along the detuning axis. Here 
= is equal to the ST + anti-crossing for every magnetic 
field, d, To convert the £* scale of the detuning axis into the s 
scale, where e = equals the 5(1, l)-5(0, 2) anti-crossing, the 
spectrum is sheared by the Zeeman-energy using the electron 
g-factor that is independently determined as explained in the 
main article. At low magnetic fields this transformation of 
scales is wrong and gives rise to a curvature of the spectrum, 
which is a function of the tunnel coupling. 



convert the £*-scale to the e-scale found in literature, for 
which £ = is defined by the 5(1,1) to 5(0,2) anti- 
crossing. To do so, we additionally shift all data points 
at a given B by \g\ns\B\ towards positive detuning (Fig. 
Sid). However, since £st+ — \9^bB\ holds true only for 
\g\isB\ ^> t C: where t c is the interdot tunnel coupling, the 
detuning axis conversion fails for low magnetic fields. As 
a result the PAT resonance-lines bend towards positive e 
for B — » T in Fig. Sid and in the spectra shown in the 
main article (Fig. 2a, b). 
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Fig. S 2: Triplet resonance in the PAT spectra, a, 

Energy eigenstates of the two-electron spin states plotted 
with the same color code as in Fig. 2c of the main arti- 
cle. Microwaves are resonant to the T + (l, 1)/T°(1, 1) transi- 
tion. Spontaneous relaxation from T°(l, 1) to the metastable 
5(0, 2) enables the detection via AIqpc- b, The PAT spec- 
trum as measured with a microwave frequency of 7.5 GHz 
exhibits a nearly horizontal blue feature as visible in Fig. 2f. 
Also at this microwave frequency the line starts from the mag- 
netic field, at which we expect the electron spin resonance 
condition to be fulfilled (black dashed line) hv — g\i(B + bo). 
c, PAT spectrum recorded at 13 GHz with corresponding line 
indicating electron spin resonance. 



Appendix C: Triplet spin resonance 

In Fig. 2f of the main article, a PAT feature is observed 
that is due to a transition from the T + (l, 1) ground state 
to the T°(l, 1) excited state. The T°(l, 1) state can relax 
via spontaneous phonon emission to the singlet bonding 
state, a superposition of 5(1, 1) and 5(0, 2). This state is 
metastable, since the spontaneous phonon relaxation is 
suppressed by the small energy difference to the T + (l, 1) 
ground state (see Fig. S2a), which makes the transi- 
tion detectable by AIqpc- A peculiarity of the T + (l, 1) 
to T°(l,l) resonance is its slope in the PAT spectrum, 
which might be a result from a gradient magnetic field 
along the magnetic field direction as discussed in the 
main article. Note that in the same detuning range, 
we observe also direct PAT transitions from the T + (l, 1) 
ground state to the singlet anti-bonding state and sin- 
glet bonding state at lower and higher magnetic fields, 
respectively. 

Here we investigate the position of the T + (l,l) to 
T°(l, 1) resonance, as a function of the microwave fre- 
quency v. The Figs. S2b and S2c show raw PAT spectra 



(without any post-processing step applied as explained 
above) recorded with v = 7.5 GHz and v — 13 GHz, re- 
spectively. The dashed lines mark the magnetic field, 
at which the electron spin resonance condition hv = 
gHB(B + bo) is fulfilled. Here we use the absolute elec- 
tronic g- factor \g\ = 0.382 and longitudinal magnetic field 
offset bo = 109 mT as determined by the Figs. 3a and 
3b of the main article. Alternatively, we might use the 
T + (l,l) to T°(l,l) resonance feature to determine \g\. 
If we use the center magnetic field of this feature as the 
resonant field, we calculate \g\ = 0.384 ± 0.005 from the 
microwave frequency dependence in good agreement with 
the \g\ = 0.382 found in Fig. 3a in the main article. 



Appendix D: The Am — PAT transition 

In the main article, we discuss whether the red (Am = 
0, An = 1) PAT resonance is dominantly due to a transi- 
tion from the 5(0, 2) ground state to the 5(1, 1) state or 
to the T°(l, 1) state. It is difficult to spectroscopically 
resolve these transitions, since they differ only by the ex- 
change energy J(e). In Fig 3a of the main article, we 
use As ± (B), the difference in detuning between the red 
Am = ±1 and the Am = lines, to assign the PAT res- 
onance. Here, we support this argument by calculating 
the expected Ae ± (B) functions for both extreme scenar- 
ios: a pure 5(0, 2)/5(l, 1) and a pure 5(0, 2)/T°(l, 1) 
transition. For the calculation, we use the determined 
\g\ = 0.382, t c = 8.7 /ieV and b = 109 mT values. 
The result is displayed in Fig. S3a for a microwave 
frequency v — 20 GHz. Obviously, the scenario of a 
pure singlet transition results in Ae^{B) > Ae~(B), 
whereas As+(B) = Ae~(B) is found for a transition to 
the triplet state. In the experiment, we clearly observe 
Ae^(B) > As~ (B) and therefore the Am = transition 
is dominantly a singlet transition. 

In both scenarios, the Ae~ (B) increases non-linearly 
towards high £?, since the resonance Am = —1 be- 
comes sensitive to the curved singlet anti-crossing at 
E z « hv — t c , where E z is the Zeeman energy. This curva- 
ture is not observed in the experiment as the Am = — 1 
transition fades out at 5 ^ 1.5 T. In both scenarios, 
the Ae+(B) functions are linear, which holds true only 
if E z <C t c — hv. The choice of a high v and an ap- 
propriate £?-range allows to extract \g\ from Ae^{B) by 
a simple linear least-squares fit as demonstrated in the 
main article. 

As a final step, we analyze the splitting of the Ae^{B) 
and Ae~(B) function quantitatively. Their offsets de- 
pend upon the exchange energy J(e,t c ), which we can ex- 
perimentally vary by t c or indirectly by v. For hv — >> 2t c , 
the detuning position of the Am = resonance be- 
comes strongly affected by J. As shown in the in- 
set of Fig. S3b, the 5(0, 2)/5(l, 1) transition (green 
arrow) shifts more towards negative detuning e than 
the potential 5(0, 2)/T°(l, 1) transition (violet arrow). 
We measure Ae+(B) for various v and determine the 
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Fig. S 3: Analysis of the Am = PAT resonance. 

a, Ae + and Ae~ calculated for the parameters \g\ = 0.382, 
t c = 8.7 /zeV, b = 109 mT, v = 20 GHz and two possible 
scenarios: The Am = line is resonant between the 5(0, 2) 
ground state 5(1, 1) (red and orange line) and T°(l, 1) (green 
and blue line). Zoom-in for low-magnetic fields in the inset. b, 
The detuning offset Asq between the transition Am = +1 for 
B — »> mT (red arrow in the inset) and the Am = tran- 
sitions from 5(0,2) to 5(1,1) (green arrow) and to T°(l,l) 
(violet arrow) , respectively, is plotted as a function of the mi- 
crowave frequency v. The error bars are determined from the 
linear extrapolation B — » mT. Shown are offsets measured 
at t c = 8.7 /xeV (filled circles) and at t c ~ 4 fieV (open rect- 
angles). The green and violet lines are least-squares fits to 
the measured offsets assuming a pure 5(0,2) to 5(1,1) and 
to T°(l, 1) transition, respectively. 



we determine a reasonable fit by assuming the Am = 
PAT resonance to be purely singlet (green line). The fit 
with a potential 5(0, 2)/T°(l, 1) transition (violet line) 
fails at v = 5 GHz. The only fit parameter used here is 
the remanence of the /imagnet, which was found to be 
w 70 mT for the fit function assuming a 5(0,2)/5(l, 1) 
transition, and « 140 mT assuming a 5(0, 2)/T°(l, 1) 
transition. The latter is very unlikely, since only a max- 
imum magnetization of bo = 109 mT was found at an 
external magnetic field of 2 T. 

As a final check, we take Ae^ values into account, 
which were determined when the DD was tuned to a 
smaller t c = 4 fieV (open squares). These data points 
cannot be fitted by the fit function that assumes a 
5(0,2)/T°(l, 1) transition at all, since the Asq /gfiB val- 
ues observed are already smaller than the remanence of 
« 140 mT, which would stay valid for the altered tun- 
nel coupling. Obviously, this leads to a contradiction, 
since all Ae~Q would become negative after subtracting 
the remanence. Only the assumption of a purely singlet 
Am = PAT transition in combination with the smaller 
remanence of ~ 70 mT, as fitted above, allows reasonable 
fitting. Our conclusion from the main article is therefore 
further supported. 

Appendix E: Spin flip-tunneling mechanism 



A^q" = Ae+(B — »• 0) by a linear fit at a sufficiently high 
magnetic field range. This procedure turned out to be im- 
practical with Aeq(B), because the Am = — 1 line fades 
out at high magnetic fields. Note that the linear extrapo- 
lation of the 5(0, 2)/T + (l, 1) transition to zero field (red 
arrow) exhibits a different detuning position than the 
5(0, 2)/T°(l, 1) transition (violet arrow), because the lin- 
ear extrapolation follows the dashed blue line in the inset 
of Fig. S3b, i.e. the linear extrapolated detuning posi- 
tion is not affected by the hybridization of the singlets. 
The detuning position of the 5(0, 2)/T°(l, 1) PAT tran- 
sition, however, is affected by the singlet hybridization, 
since it lowers the energy of the initial state 5(0,2). In 
summary, the value of As^ is always larger then zero, 
but also depends upon the nature of the Am = PAT 
resonance. 

In the inset of Fig. S3b, Ae+ for the 5(0, 2)/5(l, 1) 
transition is drawn. Obviously, Ae^ is considerably 
smaller, if the 5(0, 2)/T°(l, 1) PAT resonance dominates 
over the 5(0,2)/5(l, 1) resonance. The analysis of Asq 
is complicated by the remanence of the /imagnet, which 
is not exactly known, but should be smaller than the fully 
magnetized field of bo = 109 mT. Due to the remanence, 
the linear extrapolation of the 5(0, 2)/T + (l, 1) transition 
towards zero external magnetic field, leaves an additional 
offset on As~q . This offset, however, is independent from 
v and t c . In Fig. S3b, the extrapolated Asq values are 
plotted as a function of the microwave frequency v for two 
tunnel couplings (filled circles and open squares). Fit- 
ting the filled circles with the well-known t c = 8.7 /ieV, 



As noticed in the main text, there is a direct matrix 
element between the 5(0,2) and the (1, 1) triplet states. 
This can occur due to the nuclear spins in the barrier 
between the dots and due to spin-orbit (SO) interaction. 
We look at a toy model to examine the relative impor- 
tance of these two processes. Specifically, we consider the 
hopping matrix element for a single electron with spin 5 
between two wavefunctions associated with an electron 
on the left (\L)) and on the right (\R)) via the perturba- 
tion: 

V = j! [-a(S~ zP y - SyP~ Z ) + P( P ySy + Pz S Z )} + 

+Avo^8{r-rj)T j -5 (El) 

3 

where we have absorbed the Rashba (a) and Dresselhaus 
(/?) terms into a single spin-orbit interaction with a char- 
acteristic spin-orbit length Xso ~ 10 /im. We recall that 
m* is the effective electron mass, A w 100/ieV, vo is the 
unit cell volume and P is the nuclear spin at fj. 

We now wish to estimate the spin-flip tunneling for the 
single electron case, given by averaging over the orbital 
dipole: 

(R\ V \L) (E2) 

This can be evaluated explicitly for \L(R)) = exp(— (z ± 
a/2) 2 /4cr 2 )/(27rcr 2 ) 1 / 4 0(x, y) where z is the inter-dot axis 
(at an angle with the axis z from the spin-orbit interac- 
tion in Eq. IE1|) and <fr(x, y) is the transverse-longitudinal 
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wavefunction. As tunneling occurs only along the z-axis, 
matrix elements with p y are zero. We find two tunneling 
matrix elements: 



tso 



h 2 



particle wavefunction, cr, is related to the orbital energy 
scale of a single dot by A « — a2 . Thus, the relative 
strength of the two tunneling terms (including spin flip) 
is 



(E3) 



m*<j 2 4A50 

tnuc = #M£^nuc,f • 5 (E4) 
ft = -cos(0)[(a - (})cos(0) + (a + /3)sin(0)]£ - 
sin(0)[(/3 - a) sin(0) + (a + f3) cos(0)]£ (E5) 

9»BB nuCii = ^o^l^( r i)ll^( r i)l^' ( E6 ) 
i 

We remark that the rms value for B nuc ^ is given by 



\tso\ 

\t"nuc 



3\n x B e 



A/VN 



2|B fi - 



(E8) 



where TV is the number of spins in a single quantum dot. 

We explore briefly how this ratio varies with dot size a 
and spacing a. Specifically, A\^N ex a -1 , so a larger dot 
reduces the strength of spin-orbit tunneling compared to 
hyperfine-assisted tunneling. On the other hand, increas- 
ing the distance a increases the relative strength of spin- 
orbit tunneling to hyperfine-assisted tunneling. Setting 
gHB\j\B nnc f\ 2 = A /vg V \^j l \ 2 \^ r \ 2 I(I + 1) in A = 100 /xeV, N = 4 x 10 6 , a = 75 nm, A = 1000 /xeV 

y and A50 = 10 /im, we calculate the ratio of the matrix 



elements 



^>§5>lI 4 ( E7 ) 



[*soJ 

I tnuc I 



to be ~ 60. We remark that in external 



That is, it is the rms value for a single dot, A/vN, mul- 
tiplied by exp(— a 2 /8a 2 ). Also, the size of the single- 



field parallel to ft prevents any SO spin-charge flips. All 
spin charge flips occur only via t nuc . 

In the bases (T"(l,l), | ±t>, I U), T+(l,l), S"(0,2)) 
the two electron Hamiltonian that ignores the small t nuc 
term has the form 



/ B z 

ABj+iAB^ 
2 

AB± +iAB^~ 



ABt 



-iABt 



ABz 



iABi 





AB^r+iAB^ 





abw 

AB^+iABi 



AB^-iAB A y 

AB^-iAB^ 
2 

B 7 



2 2 ^ z 

\ -tso,y/V% (itso,z -t c )/V2 (itso,z +t c )/V2 -tso,y/V2 



—tsO,y/V2 \ 

(-itso jZ -t c )/V2 
(-itso,z + t c )/V2 

—tS0,y 

—e 



(E9) 



where, e.g., the Larmor precession frequency of an elec- 
tron spin in the left dot is B z — ABK 



Appendix F: Simulations of the PAT spectra - 
relaxation 

The simulated spectra in Figs. 2c, d of the main arti- 
cle include the effect of the phonon-mediated relaxation. 
In addition to the explanations of the simulations in ap- 
pendix A, we continue here on the coupling of the electron 
spin to the phonon bath. We thereby neglect deformation 
phonons as the energy scales examined in the experiment 
(7-22 GHz) are much smaller than the characteristic fre- 
quency scale of a phonon on the length scale of the dot 
Cph/ldot ~ 60 — 120 GHz. To determine the coupling, 
we take as an ansatz for the electronic wavefunctions the 
Fock-Darwin states, given by Gaussians, and calculate 
the coupling after orthogonalizing the states with the per- 
turbation V ph = Ee,kf( k z)^/^^f3 e (a ejk - al_ k ) 



[27j . where f(k z ) ~ 1 for the energy scales we are work- 
ing with. We then use Fermi's golden rule to calculate 
excitation and relaxation from thermal and spontaneous 
emission of phonons. Finally, we numerically solve the 
superoperator for the steady state and compare the ex- 
pectation value of 15(0, 2))(5(0, 2)| with and without the 
excitation mimicking the effect of the lock-in detec- 
tion. 



Appendix G: Measurement of spontaneous 
relaxation 

In the main article, we state that the spontaneous re- 
laxation rate from T + (l, 1) to 5(0, 2) is found to be mag- 
netic field independent. Relevant for the PAT process 
is the spontaneous relaxation at a fixed energy differ- 
ence between the T + (l, 1) and 5(0, 2) state, which is set 
by the photon energy. Thus, when changing the mag- 
netic field, the detuning e has to be changed accord- 
ingly. The measurement of the spontaneous relaxation 
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Fig. S 4: Spontaneous relaxation rate. Normalized QPC 
current averaged over r = 5/xs immediately after full mixing 
at the 5T + anticrossing for various magnetic fields at zero mi- 
crowave power. The spontaneous spin relaxation after mixing 
is measured at a distance P £ to the 5T + mixing point (see 
scheme in Fig. 2c). The dependence on the averaged time in- 
terval r is displayed in the inset for B = 1.5 T and P e = 2 mV 
together with a least-squares fit (red solid line). 



rate is done as follows: Starting from 5(0,2), we popu- 
late the T + (l,l) by 50 % via a 200 ns detuning pulse 
(28| with amplitude P £ in the absence of microwaves, 
and monitor the decay back to 5(0,2). The relaxation 
rate T s can be extracted from the time averaged lock- 
in signal AIqpc(t) = ^P- J Q r exp (— tT s ) dt, where r is 
the time spent in Pauli blockade between the pulses. 
AIo = AJqpc(O), which is independent from B, is ex- 
tracted from the fit in the inset of Fig. S4. In order to 
cover various values of the detuning and the magnetic 
field, we next fix r = 5 /is and record AIqpc as a func- 
tion of P £ for three different magnetic fields. We observe 
that AIqpc (5 /is) and thus T s are essentially indepen- 
dent of B (Fig. S4). This holds true for all P £ and hence 
for all T+(l,l)-5(0,2) energy splittings. Note that re- 
gardless of B, this energy splitting is given by P £ alone. 
This reflects exactly the situation in the PAT experiment, 
for which the microwave frequency alone sets the energy 
splitting and thus also the required phonon energy for 
the spontaneous relaxation process. 



PAT sidebands at larger detuning appear when E is 
increased as expected for PAT. Oscillation of the PAT 
amplitude as a function of E is hardly visible, because 
we cannot reach sufficiently high E and because the PAT 
lines become broadened as a function of power. One side- 
band emerges at high £* w 0.5 mV already at low E. This 
resonance stays at a constant £* for all B and is there- 
fore a transition from 5(0,2) to T+(l, 1) (Am = 1). As 
B is increased, the other transitions move towards lower 
e* while keeping the distance along e*. They are due to 
PAT transition with Am = 0. At B w 2 T, the Am 1 
PAT transition overlaps with the 2-photon Am = PAT 
transition. This is expected, since v = 11 GHz equals 
the Zeeman energy at this magnetic field in our double 
dot. 
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Appendix H: Power dependence 



In Fig. 4a of the main article, the power dependence of 
the PAT sidebands is shown for two magnetic fields. In 
Fig. S5a-e, v = 11 GHz-spectra measured with a series of 
magnetic field values are plotted, to ease keeping track of 
the resonances as they shift in detuning with the external 
magnetic field. Here, the spectra are plotted such that 
s* = is the 5T + anti-crossing for all B (compare Fig. 
Sic). The reference peaks due to pulsing to the 5T + 
anti-crossing are all aligned at £* = P £ = 1.5 mV, which 
is the pulse amplitude. The amplitude of the microwaves 
E is estimated from the attenuation of the high-frequency 
circuit at room temperature. 
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Fig. S 5: Power-dependence of the PAT spectra for different magnetic fields, a-e, Microwave induced change of the 
QPC current as a function of the double dot detuning e and the microwave amplitude E at the sample for 11 GHz. From a 
to e, the external magnetic field is increased by 500 mT. The detuning axis is calibrated so that £* = mV is equal to the 
detuning position at the S(l, 1)-T + (1, 1) anti-crossing (see Fig. lc for the calibration method). The pulse amplitude of the 
reference line is P £ — 1.5 mV. 
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